Onset of thermally driven self-motion of a current filament in a bistable 

semiconductor structure 
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We perform an analytical investigation of the bifurcation from static to traveling current density 
filaments in a bistable semiconductor structure with S-shaped current-voltage characteristic. Joule 
self-heating of a semiconductor structure and the effect of temperature on electron transport are 
consistently taken into account in the framework of a generic reaction-diffusion model with global 
coupling. It is shown that the self-heating is capable to induce translation instability which leads 
to spontaneous onset of lateral self-motion of the filament along the structure. This may occur in a 
wide class of semiconductor structures whose bistability is caused by impact ionization due to the 
negative effect of temperature on the impact ionization rate. The increment of the translation mode 
and the instability threshold are determined analytically. 

PACS numbers: 05.70.Ln,72.20.Ht,85.30.-z 



I. INTRODUCTION 



Current density filament in a bistable semiconductor 
structure with S-shaped current-voltage characteristic is 
a high current density domain embedded in a low current 
density environment Typically 
filaments appear spontaneously due to spatial instabil- 
ity of the uniform current flow when the average current 
density corresponds to the falling branch of the current 
voltage-characteristic 0- A single filament survives af- 
ter the transient (the-winner-takes-allprindple) whereas 
multifilamentary states are unstable % ^ ^ "ff, 
A static current filament may undergo further bifurca- 
tions related to teniporal and spatial instabilities 0, 0, 
EllIllIllIllIllilllllllSElEl. Spontaneous 
onset of the lateral movement is an example of such a 
bifurcation which has been recently studied both exper- 
imentally and theoretically [H Ellll H il This 
effect is of significant practical importance because the 
filament motion delocalizes the Joule heating of a semi- 
conductor structure operated in the current filamentation 
regime, and thus dramatically reduces the probability of 
thermal failure [Hill 120*, 21]. Travehng current density 
filaments bear deep similarities to other traveling spots 
in active media 0, El El IH HHlil, in particular, 
to traveling wave segments in feed-back controlled light 
sensitive Belousov-Zabotinsky reaction 29, 30]. 

The latest exerimental observations of moving fila- 
ments 0, 0, El 1^ have been made using a recently 
developed interferometric mapping technique Ex- 
periments show that in a device with a stripe-like geom- 
etry the filament travels in a soliton-like manner with a 
constant velocity ~ 1 m/s and reflects from the device 
boundaries. Both measurements and full-scale numerical 
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simulations of the transport processes suggest that the fil- 
ament motion is thermally driven and related to the Joule 
self- heating 0, ITM I21I ]. Thermally driven self-motion 
is possible due to the negative effect of temperature on 
the impact ionization rate |3lL Is^ and may occur in a 
wide class of semiconductor structures whose bistability 
is caused by impact ionization 0, E3 ■ The temperature 
acts as an inhibitor which provokes translation instability 
of a steady filament and causes its self-sustained motion 
towards a cooler region [isl E3 ■ 

In this article we present an analytical investigation of 
the translation instability which occurs at the onset of 
thermally driven filament motion. Using a recently sug- 
gested three-component reaction-diffusion model for cur- 
rent filamentation in presence of the Joule self-heating 
|22| ]. we perform a stability analysis of a static filament. 
We show that the translation invariance of a static fila- 
ment in a large structure breaks in presence of the self- 
heating. This occurs with increase of the thermal relax- 
ation time which depends on the heat capacity of the 
semiconductor structure and the efficiency of the exter- 
nal cooling. The instability threshold decreases when the 
effect of temperature on the vertical transport becomes 
stronger and increases with heat diffusion coefficient. 



II. THE MODEL 

Current filamentation in bistable semiconductors can 
be described by a two-component activator-inhibitor 
model which consists of a nonlinear reaction-diffusion 
equation, accounting for the internal activator dynamics 
in a bistable semiconductor, and an integro-differential 
equation, accounting for the applied voltage. This model, 
originally developed for bulk semiconductors with so 
called overheating instability '3, ^ and layered semicon- 
ductor structures 5, 6, 7] , proved to be efficient for a wide 
class of bistable semiconductors with various mechanims 
of bistability (see Ref. E^ and references therein). We 



2 



refer to Refs. |3J, |35|, |3a for comprehensive general for- 
mulation of this model. 

To decribe the thermally driven motion of a current 
density filament, the effect of temperature T on the 
cathode-anode transport and the heat dynamics in the 
structure should be additionally taken into account. The 
respective extension of the conventional two-component 
model (3^ |35|, ,36j has been recently suggested in Ref.22; 

da 

— =V^{Da{a)V^a) + f{a,u,T), (1) 
dT 

TT ^ = 4 A±T + {Juh + Te^t - T) , (2) 

Tu— = Uq — u — R / J [a, u) dxdy, t„ = RC. (3) 
dt Js 

Here a acts as an activator, u and T as inhibitors. 

Eas. (|ll3|) are the same as in the conventional two- 
component model 34, 33]. The variable a{x, y, t) charac- 
terizes the internal state of the device, u{t) is the voltage 
over the device. The local kinetic function / and the 
J(a, u) dependence contain the information about the 
vertical electron transport in the cathode-anode direc- 
tion and together determine the S-shaped local current- 
voltage characteristic (Fig.l). The function f{a,u) has 
three zeros corresponding to the off, intermediate and 
on states, respectively, when the voltage u is within the 
bistability range [u^, Uth] (the upper insert in Fig.l). The 
diffusion coefficient Da (a) characterizes lateral (along the 
(x, y) plane) coupling in the spatially extended semicon- 
ductor structure. Eq.Q is a Kirchhoff's equation for 
the external circuit. R is the load resitance switched in 
series with the semiconductor structure, Uq is the total 
applied voltage, C is the effective capacitance of the sam- 
ple and the circuit, 5" is the device area of the structure 
cross-section. In the following we consider the current- 
controlled regime. 

Eq.|21) describes the heat dynamics in the structure. 
Similar to the variable a, the temperature T(x,y,t) de- 
pends only on the lateral coordinates x and y. Text is 
the temperature of the external cooling reservoir, 7 is 
heat transfer coefficient, the termal relaxation time tt 
and diffusion length £t are given by 



Tt = cpw/j, 



Kw/"f, 



(4) 



where c, p and k are specific heat, density and heat con- 
ductivity of the semiconductor material, respectively, and 
w is the device thickness in the cathode-anode direction. 

The master equation is coupled to the heat equa- 
tion via the temperature dependence of the local ki- 
netic function f{a,u,T). Since the heating suppresses 
impact ionization, the effect of temperature on the ver- 
tical transport is negative and drf < 37]. The direct 
dependence of the current density J on T is neglected. 

The cha racteristic size of a filament wall is £a = 
\/Da/daf 34, 35]. Parameters £a and £t are diffusion 
lengths of the activator and inhibitor, respectively. For 
simplicity, in the following we take Da{a) = const. Next, 



we take into account only one lateral dimension x, assum- 
ing that due to the stripe geometry of the semiconductor 
structure {L^ ^ £a ^ Ly) the current density distribu- 
tion is uniform along the y direction (the lower insert in 
Fig.l). We also assume that the filament is located far 
from the edges of the stripe x — ±Lx/2. Current density 
and temperature profiles in a static filament are shown 
in the upper panel of Fig. 2. For £a ^ £t the tempera- 
ture profile essentially follows the current density profile 
(Fig. 2, curve 1 in the upper panel). For £a <C £t the 
characteristic size of the hot area is much larger than the 
filament size (Fig. 2, curve 2 in the upper panel). 



III. MECHANISM OF TRANSLATION 
INSTABILITY 

Antisymmetric (with respect to the center of the fila- 
ment) fluctuations of the current density do not change 
the total current and therefore are not suppressed by the 
external circuit in the current-controlled regime. By in- 
creasing the current denisity on one side and decreasing 
on another such fluctuation effectively leads to a small 
shift of a filament, moving one of the filament edges to a 
cooler region and another to a hotter region. Due to the 
heat inertia of the semiconductor structure the tempera- 
ture does not follow the current density instantly. Hence 
favorable conditions for the impact ionization (produc- 
tion of the activator a in terms of Eq.(Q) appear at the 
one of the filament edges, whereas at another edge the 
impact ionization is suppressed. This leads to further 
shift of a filament and, for a sufficiently long delay of the 
temperature response, to a self-sustained filament mo- 
tion. Experiments [TsL I2H show that in a high quality 
structure filaments start to move to the left and to the 
right with equal probability. This suggests that motion 
is triggered by nonequilibrium fiuctuations of the current 
density. 



IV. SETTING UP THE STABILITY PROBLEM 
A. Linearization and the eigenvalue problem 

Let us consider a stationary solution ao(a;), Tq{x), uq 
of the equations H1I2I3() : 



Daa'o + /(ao, uq, To) 0, 



rpll 



J(ao,-»o)Mo 
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rn = o, 



U^-Uq- RLy {J{aQ,uo)) = 0. 



(5) 
(6) 
(7) 



Here the prime (...)' denotes the derivative with respect 
to X and angular brackets (...) denote integration over 
the interval [—Lx/2, Lx/2]. 

Linearization in the vicinity of the stationary solution 
with respect to the fiuctiations 

a{x,t) — ao{x) = 5a(x) exp(Ct), 
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u(t) — uq = Su exp(Ct), 
T{x,t)~To{x) = ST{x) exp{(t) 

leads to the eigenvalue problem 

CSa = HaSa + dufSu + drfST, (8) 

^XT^ B , ^odaJ . , J{ao,uo) +uoduJ , . 

ttQoI = Ht ^ oa-i ou, (9) 

7 7 

T^(Su = -{l + RLy{duJ))Su- RLyidaJSa). (10) 
Here 

and all derivatives are computed at the steady state 
aoix), To(x), uo- 

B. The reference case T — Text • 

Let us briefly summarize the results of the stability 
analysis in the reference case of constant temperature 
T = Text 1^ ■ In this case Eq. is reduced to 

QSa^ Ha5a + duf 5u. (12) 

The eigenfunctions 'Si, i — 1,2,.. of the self-adjoint 
operator Ha correspond to eigenmodes of the station- 
ary filament in the voltage-controlled regime u — const 
(Fig. 2, lower panel). These eigenmodes are orthogonal 
and obey the oscillation theorem which says that '^k has 
{k — 1) nodes inside the interval [—Lx/2^Lx/1]. Due to 
the translation invariance for an infinitely large sample 
the spectrum of Ha contains a neutral (so called Gold- 
stone) mode 

*G - flo (13) 

with \q — 0. This mode corresponds to the translation 
of the filament along the x axis. For a solitary filament 
ao(a:) has a single maximum. Hence has one node 
and is identified as '^2 ^38]. Since A2 = and Ai > A2, 
the first eigenvalue Ai is positive. The first eigenmode 
^'i corresponds to spreading or shrinking of the current 
filament. 

Since ao is symmetric with respect to the center of 
the filament, the "potential term" daf in the operator 
Ha is also symmetric. Hence ^'i and \E'2 are symmetric 
and antisymmetric, respectively. The increase of ^'i is 
prevented by the global constraint Q which effectively 
fixes the total current when the external resistance R 
is sufficiently large 'sB']. In contrast, for 5a = '^2 the 
last term in the Ea. (|10|l vanishes due to the asymmetry 
of ^2- Hence in this case 5a and 5u are not coupled, 
and we can take 5u = Q. Physically, it happens because 
^'2 does not change the total current and therefore the 
voltage u also does not change. Thus for the current- 
controlled regime the dynamics of 5*2 is the same as in 



the voltage-controlled regime, and C = A2 = Ulllll. 
We conclude that for the two-component model filaments 
have neutral stability with respect to translation when 
the lateral dimension of the semiconductor structure is 
large (L^; > ^a)- 



C. Translation instability 

We shall focus on the second eigenmode '^2 with A2 = 
and investigate whether for the extended model H1I2I3|) 
the respective increment C departs from zero due to the 
coupling between the master equation and the heat 
equation As we argued above, in linear approxima- 
tion the antisymmetric translation mode does not inter- 
act with external circuit and hence 5u = 0. The set of 
equations (|8I9I1U|) reduces to 

C5a - Ha5a + dTf5TG, (14) 
ttC5T = Ht5T + {uadaJ h) 5a. (15) 

For the model (|1I2I3|I the translation Goldstone mode 
{5aG,5TQ) has two components which are given by 

5aG = ftp, 5Tg = Tq. (16) 

By taking the space derivative of Eas. l5l6|l we find that 
{5aG,5TG) indeed satisfies Eas. H14ll5ll for C = 0: 

Ha5aG + drf 5Tg = 0, (17) 
Ht5Tg + (uodaJ/j) 5aG = 0. (18) 

V. CHARACTERISTIC EQUATION 

The increment ( of the translation mode can be found 
from Eas. H14ll5|l by means of the perturbation theory 
with respect to the isothermal case — X2 — 0. Substi- 
tuting 5a = 5aG in Ea. (|15|) wc obtain 

ttCST = Ht5T + (uodaJ/i) 5aG. (19) 

The solution 5T [C, 5aG] determines the temperature re- 
sponse to the variation 5aG of the activator distribu- 
tion and depends on both 5aG and ^. Substitution of 
5T [C, 5aG] into Eq.^ leads to 

C5aG^Ha5aG + dTf5T[C,5aG]. (20) 

Now note that according to Eas. l|18ll9|) 5Tg can be pre- 
sented as 

5Tg = <5r [C - 0, 5aG] ■ 
Therefore according to Eg. ljlTI) 

Ha5aG + drf 5T[C = 0, 5aG] = 0. (21) 
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Substracting H21|l from (|20|l . muliplying by Sac and in- 
tegrating over X, we come to the equation for C: 



c 



{drf {ST[C, Sag] - 5T[C ^ 0, Sag]) a'p) 



(22) 



In the next two sections we find ST[(^, Sag] and solve 
the equation (|22|1 for the hmiting cases ia ^ £t and 
ia « ir- 



VI. WEAK HEAT DIFFUSION 

In this section we consider the case £a ^ when lat- 
eral diffusion of the internal parameter a is much more 
efficient than the heat diffusion. In this case the temper- 
ature profile essentially follows the current density profile 
(Fig. 2, curve 1 in the upper panel). In the limiting case 
£t the solution of Eq. 1)19(1 is given by 



ST[C,Sag] 



UodaJ 



lirrC + 1) 
With this input F,q.^ yields 



(5a, 



(23) 



C 



C 



C + Tj 



A, 



7 (K)^) • ^ ' 



Two solutions C — and C ^ A — t^^ of this equation 
appear to intersect at tt — in the degenerate case 
{Xg = 0) under consideration. The relevant branch can 
be chosen by starting from Ag 7^ and taking the limit 
Xg ^ (see Appendix). This leads to the piecewise 
linear dependence: 



C-A 



C = o 



for 
for 



TT < A 1; 



(25) 



According to H25|) . for sufficiently small tt the filament 
remains neutral with respect to translation. The bifur- 
cation from static to traveling filaments occurs with in- 
crease of tt at Tt = A~^ . 

Note that in the limiting case tt — the temperature 
T instantly follows all changes of the current density. In 
this case T is a local function of the parameter a and can 
be eliminated by re-defining the local kinetic function 
according to 



/(a, u) = / a, u, T = Text + 



J(a, u) u 
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(26) 



In this way the three-component model is reduced back 
to the two-component model which with respect to the 
filament stability is equivalent to the isothermal model 
for T = Toxt- Hence for tt — > filament has neu- 
tral stability, in agreement with Ea. lf^ . In the oppo- 
site case of slow self- heating {tt — 00) the increment ( 
has the largest value C = A. Generally, the heat diffu- 
sion smoothes the temperature response ST [C, Sag] , thus 
shifting the instability threshold to larger values of tt- 



VII. STRONG HEAT DIFFUSION 

In the opposite limiting case £t ^ (a the lateral 
spreading of heat is more efficient than the lateral spread- 
ing of the current density (Fig. 2, curve 2 in the upper 
pannel). We restrict the analysis to a filament with a 
flat top (Fig. 3). Such filaments are typical for large struc- 
tures {Lx ^ la). The filament width and characteristic 
width of the filament wall are denoted as W and £f, re- 
spectively. Note that if ^ L [3ll35l| . 



In this case the translation mode 



is distinct 



from zero only within the filament walls where it has a 
characteristic value ±{aon — o,off)/£f (Fig. 3, lower panel). 
Taking into account the scale separation if <^ £tj we can 
approximate Sag as 



Sag = ^fijl2 [S{~WI2) ~ S{W/2)] , (27) 

where 5{x) is the delta-function, the center of the fila- 
ment is taken at x = and the prefactor ^ if jl provides 
normalization ((Jog)^) = 1- For Sag given by H27|) the 
solution ST [C, Sag] of Eq.lO is 



ST = 
for 

ST = 
for 

ST = 
for 
where 

M(C) 



Mir, 



1 — exp 



W 

It 



exp 



x+_(W/2)\ 
It J 



(28) 



X < -W/2, 
W 
2It 



—A 'Ht exp 



sinh 



^W/2 <x < W/2, 



Mir, 



1 — exp 



W 

It 



exp 



-iW/2) \ 

iT ) 



X > W/2, 



UodaJ 



2 47(1 + ttC)^ 



iT = 



Vl + rrC 



(29) 



Substituting ST[C, *g)] into 1(22) we obtain 



C = 9t/ 




M{C)iT{C) 



^M{Q)ij 



1 — exp 



W 

'mo. 



w 



.(30) 



Taking into account that near the bifurcation point C, is 
small, we expand the right handside of Ea. H30() over (r^C) 
and come to the equation 



ttC 
2 + ttC 

AiW) 



n^^-L^^^M^A{Wl (31) 



iT 27 



w 



Ea. (|31|l coincides with Ea. (|24|l when A is replaced by 0. 
and Tt by {tt/2). Hence the solutions have the same 
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structure as in Ea. H25|) : 



for TT < 2fl-\ (32) 
for TT > 2fl-^. 



VIII. DISCUSSION 



Onset of filament motion for weak and strong 
heat diffusion 



For £t <C ia the condition for the onset of filament 
motion is given by Eg. 1)25(1 and can be presented as 



TT > Ti 



7 



uo {dTf da J (a'^) 



(33) 



In the opposite hmiting case £t S> ia this condition is 
given by 



TT > T2 = 



Ail 



7 



tfA{W) UodTfdaJ' 



(34) 



In both cases £a ^ and £a ^ It the instabihty occurs 
when the thermal relaxation time tt is sufficiently large. 
This suggests that tt is a bifurcation parameter also in 
the intermediate case £a ^ £t- 

The relation between t^^ and r^*^ becomes transparent 
when OtJ and da J are considered as constants. (Note 
also that £f ^ £a-) In this case 



th 



A£i 



til 



laA{W) 



(35) 



For wide (VK ^ £t) and narrow {yV <C £t) filaments the 
expansion of yl(T/F) over WjiT (see Eq.(|3TJ) leads to 



-th 



-th 



£t_ 

' Ta 

2{£Tf 

£aW^ 



.th 



-th 



for 
for 



W->£t, 
W < £t. 



(36) 



Thus the threshold time t^ exceeds t'^ by a factor which 
is always larger than £T/£a- 

For £t ^ £a the instability threshold has been previ- 
ously obtained in Refl23 for the current density profile 
as shown in Fig. 3. The argument is based on the find- 
ing that the velocity v of the stationary filament motion 
is proportional to the difference of temperatures at the 
front and back filament edges AT: v = CAT 0|. This 
temperature difference, in turn, depends on the v. Hence 
the velocity of stationary motion can be found from the 
equation v — CAT{v). Since for v = the temperature 
profile is symmetric, AT(0) = 0. Therefore w = is al- 
ways a solution of this equation. However, this solution 
is unstable if AT{v) increases faster than v/C with in- 
crease of V. In this case when v deviates from v = 0, the 
"overproduction" of the the temperature difference AT 



leads to further increase of v. The respective condition 
is given by 



Jotf)uo 



7 (K)2) 



OtI da. (37) 



Here vq has the meaning of the upper limit of the filament 
velocity and vt = £t/tt is the thermal velocity. Jon , 
Joff and Ooff , flon are the current denisties and the values 
of the variable a correspoding to the uniform on and off 
states at u — uq, respectively (Fig. 3). ^1(1^) is exactly 
the same as in 

Eq.lISni- The criterium ^ can be re- 
written as 



TT > T2 , 



;sth 



(38) 



A{W) vo ■ 

Taking into account that for a wide filament 

((a[,)2) «2(ao„-aoff)V^/ 
and considering the function Ot/ as a constant, we find 

5t/ "0 Jon — Joff 



(39) 



27 Oo 



Ooff 



Hence 



crth 
T^2 



floff 



£fA{W) UodTfJon-Jofi 



(40) 



(41) 



Comparing Ea. H34|l and Ea. (|41|l . we see that T2^ coin- 
cides with T2^ (see Eg . (1341) if we approximate the partial 
derivative da J by the finite difference 

daJ ^ (Jon - Joff)/ (aon " floff)- 



B. Traveling filament 

Remarkably, the difference between the regimes of 
weak and strong heat diffusion, which is well-pronounced 
near the bifurcation from static to traveling filament, 
vanishes for fast self-sustained filament motion with ve- 
locity V ^ Vt- In this regime the static thermal diffu- 
sion length £t = \/Dttt should be replaced by Ct = 
^DTTf^ where r/ — W/v is the time of the filament 
passage over its width W and Dt = n/ cp is thermal dif- 
fusivity. For Ct ^ W the heat diffusion is negligible. 
This condition yields 



Dt/W = vt£t/W. 



(42) 



For typical filament width W ^ IQ iim the condition 142() 
yields v 3> 10'^ cm/s for Si and v ^ 250 cm/s for GaAs 
{Dt — 0.92 and 0.25 cm^/s, respectively). 

If the condition Ct ^ is met togeter with the con- 
dition Tf ^ Tt, the heating is not only local but also 
adiabatic. (Note that for £t ^ W these two conditions 
are equivalent.) In this case the rate of temperature in- 
crease in each point is proportional to the local power 
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dissipation. The velocit y o f such fast narrow filament 
has been obtained in Refl23: 



V W y/vo W/tt- 



(43) 



Substituting the definitions tt and vq (see Eas. (l4l40ll V 
we find that the filament velocity is proportional to the 
square root of the total Joule power which is dissipated 
in the filament: 



K ■ P 

v = J , P 

V cp 



(Jon - Jos)Wuq JonWuQ 



IV 



W 



(44) 



where cp is specific heat per unit volume and the coeffi- 
cient K is given by 



K 



1 



drfda - 
Qoff 2(a, 



Ifdrf 



floff j 



(45) 



Alternative elementary derivation of the square root de- 
pendence is presented in Ref,21. 

Apparently, the square root dependences H43I44|) re- 
flect the experimental situation which is favorable for 
observation of traveling filaments because large veloc- 
ity corresponds to low instability threshold. It has 
been reported j23| that the results of experimental mea- 
surements of the filament velocity in a Si device {v ~ 
3 • lO^'cm/s, W ~ 10 fim) indeed obey the square root 
dependence. 



C. The type of bifurcation 



The piecewise linear dependence H25I32|I of the incre- 
ment ( on the bifurcation parameter t^^ as well as the 
square root behaviour of the filament velocity near the 
bifurcation point, predicted in Refl22L are qualitatively 
the same as obtained in Refl2^ for the three-component 
"cubic" model for a traveling spot. These features seem 
to be typical for a degenerate "drift pitchfork" bifur- 
cation associated with Xq = 0. Note that the mecha- 
nism of self-motion described here for one-dimcnsional 
spatial domains is also valid for two-dimensional do- 
mains. The difference between the model (|1I2I3II and 
three-cornponent models for traveling spots studied in 
Refs-imililliill^ has been discussed in Ref^ 



For this mechanism inhibition is associated with the in- 
ternal voltage drop across the plasma layer inside of the 
semiconductor structure. The regions of activation and 
inhibition are spatially se par ated. The effective reaction- 
diffusion model 0, 14^ l43l | for such two-layer system 
has the same structure as (|1I2I3|I . Recent numerical 
simulations demonstrate that this mechanism leads to 
filament motion during the switching-off transients in 
power diodes The applicability of this concept to 

thyristor-like — n+ —p — n~ — n~^ structures discussed 
in llol should be considered in view of the latest stud- 
ies |44l I45I |46| which reveal complexity of the vertical 
transport processes in these multilayer structures. 



IX. CONCLUSIONS 

Joule self-heating of a static current-density filament 
may lead to the translation instability and the onset of 
thermally driven self-motion. This may occur in a wide 
class of semiconductor devices whose bistability is caused 
by impact ionization because impact ionization coeffi- 
cients decrease with temperature. The eigenvalue of the 
respective translation mode, which is zero when heat- 
ing is neglected, becomes positive due to the suppressive 
effect of temperature on the electron transport in the 
cathode-anode direction. Increments of the translation 
mode are found for the limiting cases of weak (|25l) and 
strong H32() heat diffusion. In both cases the instability 
threshold (|33I34|) is controlled by the thermal relaxation 
time Tt and is directly proportional to drf- It scales as 
Irl^a with increase of the thermal diffusion length It- 
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APPENDIX: INCREMENT OF THE 
TRANSLATION MODE 



D. The alternative mechanism for self-motion of a 
current filament 

Generally, the self-motion of a current filament be- 
comes possible when a semiconductor structure, apart 
from the activator mechanism which leads to bistabil- 
ity, posseses an internal mechanism of slow inhibition. 
As we described here, such mechanim can appear due to 
the effect of Joule heating on the cathode-anode electron 
transport. Another inhibitor mechanism of purely elec- 
trical origin has been discussed in Refs.fT^. 1391 llol |4]| . 



To indentify the physical branch of the solution of 
Ea. (|24|) . we note that for Ac 7^ this equation becomes 



C 



A. 



Smooth solution of this quadratic equation 

c 



(A.l) 
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for Ag ^ tends to the piecewise linear dependence 
given by F^q.^^. Note the similarity of Ea. (|A.llA.2() 



and Eq.(9) in Ref[2E 
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FIG. 1: S-shaped current-voltage characteristic J{u) of a bistable semiconductor structure. The upper insert shows the local 
kinetic function f{a,u) for u < Uh, Uh < u < uth and u > Uth (curves 1,2 and 3, respectively), where Uh and Uth are the 
hold and the threshold voltages at the edges of the bistability range. For Uh < u < Uth the local kinetic function has three 
zeros corresponding the off, intermediate and on branches of the current-voltage characterisics, respectively. The lower insert 
sketches the rectundular semiconductor structure elongated along the x axis. 




FIG. 2: Profile ao(x) of the activator and the temperature profile To{x) in a steady filament (the upper panel). The current 
denisity profile J{ao{x),u) is quailitatively the same as ao{x). The temperature profile To{x) is sketched for weak {£t <S £a, 
curve 1) and strong {£t S> £a, curve 2) heat diffusion. Note that in our analysis Lx is much larger than the filament width. 
The lower panel shows two first eigenfunctions ^'i and ^'2 of the operator Ha and their eigenvalues Ai, A2. 'I'l and Vl/2 are 
eigenmodes of the steady filament in a voltage controlled regime and correspond to expansion(shrinking) and translation of the 
filament, respectively. In the isotheramal case A2 = and ^'2 = ^'g ~ o-'q. 
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FIG. 3: Profile ao{x) of the activator and the temperature profile To (a;) in a wide steady filament (the upper panel). The 
filament width W , the width of the filament wall if and the thermal diflfusion length It are indicated. Note that in our analysis 
Lj: W £f and £t S> if, whereas the £t/W ratio is arbitrary. As in Fig. 2, the low panel sketches two first eigenmodes "ifi 
and 5'2 of the operator Ha and their eigenvalues Ai,A2. 



